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Abstract: Supersymmetric models with t—b—r Yukawa unification at Mqut qualitatively 
predict a sparticle mass spectrum including first and second generation scalars at the 3- 
15 TeV scale, third generation scalars at the (few) TeV scale and gluinos in the sub-TeV 
range. The neutralino relic density in these models typically turns out to lie far above the 
measured dark matter abundance, prompting the suggestion that instead dark matter is 
composed of an axion/axino mixture. We explore the axion and thermal and non-thermal 
axino dark matter abundance in Yukawa- unified SUSY models. We find in this scenario 
that i). rather large values of Peccei-Quinn symmetry breaking scale fa ~ 10^^ GeV are 
favored and ii). rather large values of GUT scale scalar masses ^ 10 — 15 TeV allow for the 
re-heat temperature Tr of the universe to be Tr ~ 10^ GeV. This allows in turn a solution 
to the gravitino/Big Bang Nucleosynthesis problem while also allowing for baryogenesis 
via non-thermal leptogenesis. The large scalar masses for Yukawa-unified models are also 
favored by data on ?j — > 57 and Bs — > fJ''^t^~ decay. Testable consequences from this 
scenario include a variety of robust LHC signatures, a possible axion detection at axion 
search experiments, but null results from direct and indirect WIMP search experiments. 
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1. Introduction 

The celebrated unification of gauge couplings at scale Mqut — 2 x 10^^ GeV under mini- 
mal super symmetric standard model (MSSM) renormalization group evolution|j^] strongly 
suggests that nature is described by some sort of supersymmetric grand unified theory 
(SUSY GUT) model at very high energy scales. While the GUT gauge group ^[/(S)^ has 
many compelling and beautiful properties, the gauge group SO (10) has an even greater 
appeal]^, and also some indirect experimental support in terms of how well see-saw neu- 
trino mass fits into the general scheme|Q]. While both SU{5) and SO{W) SUSY GUT 
theories admit gauge coupling unification, 50(10) theories also yield matter unification, in 
that all superfields of a single Standard Model (SM) generation- plus a SM gauge singlet 
superfield A^'^ containing a right-hand neutrino state- fit neatly into the 16-dimensional 
spinorial representation V'le of 50(10). 

In the simplest 50(10) SUSY GUT models, the MSSM Higgs superfields- Hu and Ha- 
both live in the fundamental representation ^iq. In these models, the superpotential has 
the form 

/ 9 /V'leV'ie'^io + • • • (1-1) 

where the dots represent model-dependent terms which, for instance, might include further 
Higgs fields responsible for the GUT gauge symmetry breaking. The coupling / represents 
the unified Yukawa coupling of each generation: thus, just as 5C/(5) models often predict 
fb — fr unification, simple 50(10) SUSY GUT models predict the more restrictive t — b — r 
Yukawa coupling unification at the GUT scale |^]. 

It is probably fair to say that at this moment no compelling 50(10) SUSY GUT 
model yet exists. Models based in four spacetime dimensions require large, unwieldy Higgs 
representations to break the 50(10) GUT symmetry. Newer models formulated in extra 
spacetime dimensions are able to do away with the large Higgs reps and break the GUT 
symmetry via compactification of the extra dimensions In our approach here, we will 
adopt a pragmatic view, assuming that the MSSM (or MSSM plus gauge singlets) is the 
correct effective theory describing physics between the weak and GUT scales, and we will 
explore some of the consequences of requiring the three third generation Yukawa couplings 
to have a high degree of unification at Mqut, as suggested by simple 50(10) SUSY GUT 
models. 

In our calculation to check third generation Yukawa coupling unification, we begin with 
the measured gauge couplings and third generation fermion masses, and use renormalization 
group methods to evolve the coupled gauge and Yukawa couplings up to the GUT scale. The 
calculation ends up being sensitive to the entire SUSY particle mass spectrum through weak 
scale threshold corrections involved in transitioning between the SM and MSSM effective 
field theories 1^]. The parameter space of the model is given by 

mi6, mio, M|,, mi/2, Aq, tan/3, sign(/i) (1.2) 

where mig is the GUT scale mass of all matter scalars, mio is the GUT scale mass of Higgs 
scalars, M^) is their D-term value, mi/2 is the (unified) GUT scale gaugino mass, Aq is the 
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unified GUT scale SSB trilinear term, tan/3 = Vu/vd is the weak scale ratio of Higgs field 
vevs, and is the superpotential Higgs bilinear term, whose magnitude-but not sign-is 
determined by the scalar potential minimization conditions. 

In practice, the two Higgs field soft breaking terms- mn^u^ and "i^^- cannot be de- 
generate at Mqux and still allow for an appropriate radiative breakdown of electroweak 
symmetry (REWSB). Effectively, rn?ii^ must be less than tn^n^ at Mqut in order to give 
rn^H^ a head start in running towards negative values at M^eafc- We parametrize the Higgs 
splitting as m?^ ^ = m^Q =p 2M|, in accord with nomenclature for L'-term splitting to 
scalar masses when a gauge symmetry undergoes a breaking which reduces the rank of 
the gauge group. A D-term splitting should apply to matter scalar SSB terms as well; in 
practice, better Yukawa unification is found when the splitting is only applied to the Higgs 
SSB terms. Such a GUT scale Higgs mass splitting might arise via GUT scale threshold 
corrections!^. 

In previous work, the above parameter space was scanned over (via random scans P, 10 1 



and also by more efficient Markov Chain Monte Carlo (MCMC) scans|ll]) to search for 



Yukawa unified solutions using the Isasugra subprogram of Isajet|12| for sparticle mass 
computations. The quantity 



j^^ max{h,hJ ) (evaluated at Q = Mgut), (1.3) 
mmiftjbjr) 



was examined, where solutions with i? ~ 1 gave apparent Yukawa coupling unification. 
For superpotential Higgs mass parameter > (as favored by {g — 2)^ measurements), 
Yukawa unified solutions with i? ~ 1 were found but only for special choices of GUT scale 
boundary conditions P, |l^, B, ITOI, 111, |ll|: 



mi6 ~ 3 — 15 TeV, ~ — 2mi6 , miQ ~ 1.2mi6 , rni/2 ^ rriiQ , tan/3 ~ 50 . (1-4) 

Models with this sort of boundary conditions were derived even earlier in the context of 
inverted scalar mass hierarchy models (IMH) which attempt to reconcile suppression of 
flavor-changing and CP- violating processes with naturalness via multi-TeV first/second 



generation and sub- TeV scale third generation scalars|l5|]. The Yukawa-unified spectral 



solutions were thus found in Refs. [|^, 11 1 to occur with the above peculiar choice of 



boundary conditions as long as miQ was in the multi-TeV regime. 



Based on the above work|10|, 11 1, the sparticle mass spectra from Yukawa-unified SUSY 



models are characterized qualitatively by the following conditions: 

• first and second generation scalars have masses in the 3 — 15 TeV regime, 

• third generation scalars, fi and have masses in the TeV to few TeV regime (owing 
to the inverted scalar mass hierarchy). 



gauginos, including the gluino, have sub- TeV masses, 

the lightest neutralino Xi is nearly pure bino with mass typically m^o ~ 50 — 80 GeV. 
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The presence of a bino-hke Xi along with multi-TeV scalars gives rise to a neutralino cold 
dark matter (CDM) relic abundance that is typically in the range Q.^oh'^ ~ 10 — 10^, i.e. 
far above the WMAP measured] 17] value 

QcDAih'^ = 0.110 ± 0.006 (1.5) 

by several orders of magnitude. 

One solution to the CDM problem in Yukawa-unified SUSY models occurs for mig ~ 
3 TeV, where one can have m^o ~ mh/2. In this case, the x? would pair-annihilate in the 
early universe through the light Higgs h resonance at a sufficient rate to obtain the desired 



relic density|ll|, and could hence be the stable lightest SUSY particle (LSP). However, 
since miQ is relatively low, Yukawa coupling unification only occurs at the R ~ 1.09 level. ^ 
Another very compelling way out of the Yukawa-unified dark matter abundance prob- 
lem occurs if one makes the additional assumption that the strong CP problem is solved 
by the Peccei-Quinn mechanism] 18 1, which leads to the presence of a light pseudoscalar 



particle, the axion a [T^, |2^. Since we are working in a supersymmetric theory, the axion 
occurs as part of an axion supermultiplet|2l|, which contains not only the axion, but a spin- 
saxion (with mass of order the weak scale), and a spin-^ axino a. The axino is i?-parity 
odd. While the saxion is expected to have a mass of order the SUSY breaking scale, the 



axino mass is very model dependent, and can lie anywhere in the keV-GeV range[22, 23]. 
The axino then can serve as the LSP instead of the lightest neutralino |2^, |25|] . 

In the case of an axino LSP, the supposed neutralino relic abundance is greatly reduced 
since x? ^ ^7 decay can occur with a lifetime in the range of ~ 10~^ — 10^ sec (depending on 
parameters). This decay time is sufficiently short that late-time neutralino decay to axino 
in the early universe should not upset successful predictions of Big Bang Nucleosynthesis 
(BBN)]p3|]. The neutralino abundance then gets reduced by the ratio rria/m^o which can 
be in the range 10~^ — 10~^. The d coming from x? decay would actually constitute warm 



dark matter as long as ~ 1 GeV[26]. However, axinos can also be produced thermally 



in the early universe, and these would consititute CDM as long as ~ 100 keV]| 



Thus, in the axino LSP case, we could have a mixed dark matter (DM) scenario]28] with 

• thermally produced cold or warm axino DM, 

• an admixture of warm axino DM arising from x? ^ 07 decays and 

• a possibly large presence of cold axion DM. 

The axino LSP scenario turns out to be even more compelling cosmologically than just 
as a means to reconcile the dark matter relic abundance with Yukawa- unified models. In 
this class of Yukawa-unified solutions with mie in the multi-TeV range, we expect from 
supergravity theory that scalar SSB terms should be directly related to the gravitino mass 
777-3/2, and so we also expect the gravitino G to lie in the multi-TeV range. The cosmological 
gravitino problem- wherein gravitinos produced thermally in the early universe suffer a 



^Another possibility, having rriA light enough that XiXi can annihilate through the A resonance, appears 
to be excluded because these cases violate limits on BF{Bs fi'^ decay 
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late-time decay, thus destroying the successful predictions of BBN- can be avoided. For 
^^^•3/2 ~ 5 TeV, the re-heat temperature Tr must be Tr ~ 10^ GeV[29|, thus creating 



tension with most viable mechanisms for baryogenesis[30|. However, for 771,3/2 ~ 5 TeV, 
the re-heat bound is much higher: Tr ~ 10^ — 10^ GeV. This range of Tr is too low for 
thermal leptogenesis (which requires Tr ~ 10^'^ GeV)|Q, but is exactly what is needed 
for baryogenesis via non-thermal leptogenesis pl| , wherein the heavy right-hand neutrino 
states are not produced thermally, but rather via inflaton decay. It was pointed out in 
Ref. |32U that this is also the exact range needed to generate a dominantly cold axino DM 
universe. Thus, the whole scenario fits together to offer a consistent cosmological picture 
of BBN, non-thermal leptogenesis and CDM composed of axions and/or axinos, and one 
solves the strong CP problem to boot [ 32 1! 



In this paper, we make a detailed study of the cosmological consequences of Yukawa- 



unified SUSY GUT models with an axino LSP. The initial study proposed in Ref. |32| 
assumed only a negligible component of axion dark matter. In this study, we now fold 
in the axion contribution to the dark matter abundance. In addition, we have updated 
the value of rut in our calculations, and refined the Yukawa coupling 1-loop beta function 
threshold effects in Isajet. We also incorporate in this study the results of our Markov 
Chain Monte Carlo (MCMC) approach to finding Yukawa-unified solutions. Whereas the 



Yukawa-unified solutions found via a random scan in Ref. |32] had miQ ~ 15 — 20 TeV, 
the more efficient MCMC scans used here are able to find many solutions for tuiq values 
as low as 3 — 10 TeV. 

The paper is organized as follows. In Sec. ||, we update our sparticle mass predictions 
for Yukawa-unified SUSY models from Isajet using 1. an improved beta- function threshold 
decoupling, 2. an updated value of the top mass nit = 172.6 GeV, and 3. Markov Chain 
Monte Carlo (MCMC) scans of the parameter space. We also exhibit plots of the 6 ^ 57 
and Bg fJ-^t-i-^ branching fractions versus m-ig and find these favor m-ig values ~ 10 
TeV. Moreover, we perform MCMC scans for Yukawa-unified solutions with an alternative 
spectrum generator, SOFTSUSY, and compare these results with those gained from Isajet. 
In Sec. m, we review elements of axion and axino dark matter cosmology, including plots 
of the axion relic abundance, axino relic abundance and neutralino lifetime. In Sec. ^ we 
discuss the gravitino/BBN problem, non-thermal leptogenesis via infiaton decay, and mixed 



axion/axino dark matter. Our main findings are located in Sec. O. Here, we calculate all 
three components of axion/axino dark matter in Yukawa- unified models for ruiQ values of 
5, 8, 10 and 15 TeV. We explore scenarios wherein axions constitute either most, or hardly 
any, of the dark matter abundance. Our results are presented in the rua vs. Tr plane. We 
find that models with a rather large value of the Peccei-Quinn breaking scale fa ~ 10^^ 
GeV are favored, as well as models with miQ ~ 771-3/2 on the high side: ~ 10 — 15 TeV. In 
these cases, the re-heat temperature of the universe can range above 10^ GeV, allowing for 
a solution to the BBN gravitino problem as well as allowing for non-thermal leptogenesis. 
Our conclusions are presented in Sec. H. In an Appendix, we list updated Yukawa-unified 
benchmark points from Is A JET 7.79. 
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2. Updated spectrum predictions for Yukawa-unified SUSY models 



2.1 Updated Isajet calculation of sparticle mass spectra 

The Isajet calculations begin by adopting the fermion mass boundary conditions that 
mf^{Mz) = 2.83 GeV, m^{Mz) = 1.7463 GeV, and mtipole) = 172.6 GeV, along with 
the measured gauge couplings; in particular we take ^(Mz) = 0.1172. Note the value 



of mt represents an update due to recent DO and CDF measurements] 33 1 over our previous 



work Ref. |11] which was done using mt = 171 GeV. Our results do not change qualitatively 



upon varying the fermion masses within their error bars[10|. We use the Isajet 7. 79 [12 1 
program to perform two loop RG evolution of gauge and Yukawa couplings and all soft 
SUSY breaking terms. For gauge and Yukawa couplings, Isajet actually uses an RGE 
approach wherein the 1-loop beta functions change whenever a SSB threshold is passed 
overp^, |35[ . The Yukawa coupling evolution depends on finite terms from 1-loop MSSM 
threshold effects]^, these are implemented at a scale Q = MsusY = ^/^^WJ^H^- 
threshold effects cause the entire calculation to depend sensitively on the sparticle mass 
spectrum, which enters the various loop corrections to ft-, ft and fr- In addition, in the 



two-loop RG running|37| of SSB terms from Mqut to the weak scale, non-mixing soft terms 



are frozen out at their own mass scale, while SSB terms that mix are frozen at the scale 



MsusY Il38| . Complete 1-loop corrections are then applied to all sparticle masses. This 



approach leads to good agreement with other publicly available sparticle spectra codes pQ] 
when the sparticle masses are all nearby in mass scale. For spectra suffering severe splitting 
(as will be the case here), the Isajet multiple decoupling approach attempts to deal with 
the fact that several mass scales may be present in and around the weak scale ||35[| . 

Two updates in the Isajet 7.79 code affect the running of Yukawa couplings. In early 
versions of Isajet, all squark contributions to RG running were decoupled at a common 
squark mass scale taken to be m^^ and all sleptons were decoupled at a common scale taken 
to be rrie^- In Isajet 7.79, the first/second and third generation squarks and sleptons 
decouple at the values of the corresponding soft SUSY breaking terms: thus, in a case 
where rriiQ = 10 TeV, first/second generation squarks and sleptons decouple around 10 
TeV while third generation squarks and sleptons decouple at a much lower scale around 
3 TeV. In addition, earlier versions of Isajet included 2-loop terms for MSSM running 
between MsusY and Mgut, but turned these off for Q < MsusY-, where the SM was 
the expected effective theory. In Isajet 7.79, the MSSM two-loop terms remain in for 
Q < MsusY, since the scale of decoupling of 2-loop terms is a 3-loop effect. The current 
version should give a better estimate of Yukawa coupling evolution in models with severe 
first/second and third generation splitting, as in Yukawa-unified models. 



2.2 Isajet/Softsusy comparison for Yukawa-unified spectra 

Other public spectrum codes follow a different approach than Isajet and perform a 1- 
step decoupling of SUSY particles at Q = MsusY- This "all-at-once" transition may 
lead to some differences in the sparticle masses in particular in the case of a widely split 
spectrum ||3^ . It is therefore interesting to compare results from different spectrum codes. 
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Figure 1: Evolution of ft, fb and /r from tiie wealc scale to tlie GUT scale for ttie point in Table 
ISAJET results are solid, while Softsusy results are dashed. The large jumps around 3 TeV from 
ISAJET correspond to the SM-MSSM threshold corrections. 



Here, we choose Softsusy2.0.18|40] for a representative comparison of the 1-step decou- 
pling with the multi-step approach implemented in Isajet7.79. 

The evolution of Yukawa couplings between the electroweak and GUT scales is illus- 
trated in Fig. |l] for a point with rn-ig = 8 TeV, which runs for both ISAJET and SOFTSUSY. 
As can be seen, Softsusy implements threshold corrections to the [gauge and] Yukawa 
couplings at scale Q = Mz, and then runs from Mz to Mqut with full MSSM RG evolu- 
tion. The SM-MSSM threshold corrections in Isajet occur at Q = MsusY, and give rise 
to the the discontinuities located around MsusY ~ 3 TeV. The steep slope for ft and f}, 
for Q < MsusY occurs mainly because the co-efficient of the QCD g'^ contribution to /^/^ 
running changes from 8 in the SM for Q < MsusY to 16/3 in the MSSM for Q > MsusY- 
The two programs obviously differ in the bottom Yukawa coupling /{,. Part of this dif- 
ference comes from somewhat different weak scale boundary conditions on the value of 



m 



DR 



(Mz)- as mentioned above, Isajet has a hard-coded value of mj^ (Mz) = 2.83 GeV; 
in Softsusy, on the other hand, m^^{Mz) is computed from m^'^{mb) with SUSY cor- 
rections added at M^.^ The larger value of /b in Softsusy causes fr to run to slightly 
higher values at Mqut than for Isajet . The parameters of the point used and resulting 
mass spectra are listed in Table ||. Softsusy gives Yukawa unification at 2%, Isajet 
gives Yukawa unification at 8%. While many masses are very similar, the values of rrig and 
especially niA differ quite a bit. 

We next scan over the parameter space Eq. (|l.2| ) using a MCMC algorithm, as detailed 
in [11 1, that searches for solutions with R as low as possible. Figure |2| shows the degree 
of Yukawa unification, R, found with the two codes as a function of mig. Points from 
Isajet 7.79 are shown in blue, while points from Softsusy 2.0.18 are shown in green. As 



^To comply with the bottom mass used by Isajet in other parts of the program, we use m*^ (mi,) = 
4.2 GeV as input in Softsusy, leading to mf ^(Mz) ~ 2.6 GeV depending on the exact parameter point. 



parameter 


SOFTSUSY 


ISAJET 


ft 


0.559 


0.555 


fb 


0.549 


0.512 


fr 


0.560 
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1.02 


1.08 
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bi 
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2653 




8064 
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116 
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120 


120 




53.4 


53.1 


niA 


681 


1955 




129.1 


127.5 



Table 1: Masses in GeV units and parameters for a Yukawa-unified point using Softsusy 2.0.18, 
ISAJET 7.79. We take mig = 8000 GeV, mw = 9760.3 GeV, Md = 2435.7 GeV, toi/2 = 67.8752 
GeV, Ao ^ -16007 GeV, tan/3 = 48.71 and mt = 172.6 GeV. Thus, mn, ^ 10350.3 GeV and 
niH^ = 9132.3 GeV 




4 5 6 7 8 9 10 
mi 6 [TeV] 

Figure 2: Degree of Yukawa unification, R, versus mig from MCMC scans over the 5*0(10) model 
parameter space. Isajet7.79 resuhs are in blue, while Softsusy 2.0.18 results are in green. 

one can see, with both programs a high degree of Yukawa unification can be found, and 
these solutions tend to prefer high mig. However, as is apparent from the density of points, 
solutions with low R are more easily found with ISAJET than with SOFTSUSY. Another 
important difference is that Softsusy generates Yukawa-unified solutions only for mig up 
to about 9.5 TeV, while Isajet generates solutions for mig well beyond 10 TeV (here we 
only show results up to 10 TeV). The reason lies in the fine-tuning of the ^ parameter and 
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Figure 3: Points with R < 1.05 from MCMC scans projected into various planes: (a) mi/2 vs. 
mi6, (b) Aq/ttiiq vs. toiq, (c) mio vs. mig, and (d) ^ vs. mie- 1sajet7.79 results are in blue, while 
SOFTSUSY 2.0.18 results are in green. 

the electroweak symmetry breaking. 

Overall, Yukawa-unified solutions obtained with either program obey the same qualita- 
tive conditions Eq. ([1.4|). This is illustrated in detail in Figs, ^(a)-(d), which show MCMC 
scan points with R < 1.05 projected into various planes. In particular, frame (a) shows that 
mi6 is in the multi-TeV range while ~ 100 GeV. From frame (b), which shows the 
Ao/miQ vs. mi6 plane, we see that both programs require Aq ~ —2miQ, although Isajet 
prefers somewhat more negative values of Aq than Softsusy does. Frame (c) illustrates 
the correlation mio ~ 1.2mi6. Finally, frame (d) shows an important difference in results 
from the two codes: in spectra obtained with Isajet the fi parameter grows with rriiQ 
and reaches values of a few TeV, while in spectra obtained with Softsusy, the value of ^ 
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Figure 4: Points with R < 1.05 from MCMC scans over the SO{10) model parameter space, in 
(a) VL^oh"^ vs. mi6 and in (b) fl^oh'^ vs. R. 1sajet7.79 resuhs are in blue, while Softsusy 2.0.18 
results are in green. 

remains around 1-2 TeV. 

In Fig. §(a), we show the value of neutralino relic density from both Is A JET and 
Softsusy versus mie for points with R < 1.05. The relic density is calculated with 
MiCROMEGAsI^, which easily interfaces with both Isajet and Softsusy. Here, we see 
Softsusy predicts Jl^o/i^ ~ l - lOOO, with some points extending down to Jl-o/i^ ~ 0.3. 
These latter points occur due to neutralino annihilation near the light Higgs h resonance, 
but do not quite reach the WMAP-measured dark matter density because the higgsino 
fraction fn of the x? is only few per mil. Isajet tends to give an even smaller higgsino 
fraction (larger /x parameter, c.f. Fig. |3|(d)) and hence an even larger neutralino relic density 
up to ^^oh'^ ~ 10000. Points with ^^oh'^ ~ 0.1 due to annihilation through Higgs can be 



found with both programs, but then the Yukawa unification is only R ~ 1.07 — 1.09[11| 
while Jh is of the order of 1%. A plot of ^^oh? versus R is shown in Fig. §(b). 

We conclude that both Isajet and Softsusy predict a large over-abundance of neu- 
tralino dark matter from Yukawa-unified models with R < 1.05. The postulation of ax- 
ion/axino dark matter allows one to reconcile the dark matter results with Yukawa unifi- 
cation (as well as solving the strong CP problem). 

2.3 Results for b —>■ and Bg —>■ fJ-^fJ-^ decay 

In this section, we present results for the branching fractions for 6 — > 57 and Bg M^A*^ 
decays. We adopt the Isajet Isatools|42] program for these calculations. The BF{b — > 57) 



decay rate is calculated in Isajet by evaluating the Wilson co-efficients for the relevant 
operators mediating 6 ^ 57 decay at a scale Q > MsusY, and then running down to Mz 



scale using a tower of eff'ective theories approach 1 44 1. At Mz, the Wilson co-efficients are 



matched to the SM ones, and run using 2-loop evolution down to scale m-b, where the 5 — > S7 
decay rate is evaluated, including complete NLO corrections! 45]. Large tan/3 effects are 




Ql • ' • ' ' ' • ' • ' • 1 

4 6 8 10 12 14 16 

(TeV) 

Figure 5: Predictions for BF{b 57) vs. mig from Yukawa-unified models generated by 
Isajet7.79 with R < 1.05 and mt = 172.6 GeV. The blue dots show results from MCMC scans 
with fixed mie. We also show the measured central value of BF{h 57) along with ±2cr error 
bands. 



accounted for by adopting the ISAJET running Yukawa couplings, which include threshold 



effects as noted above@. The Isajet SM result, BF{b 57) ~ 3.1 x 10~^, agrees well 
with a recent evaluation by Misiakp7||, which finds BF{b s-y)sM = (3.15±0.23) x 10"''. 

The branching fraction BF{h sj) has been measured by the CLEO, Belle and 
BABAR collaborations; a combined analysis finds the branching fraction to be BF{b — > 
57) = (3.55 ± 0.26) X 10""^. In Fig. |, we plot BF{h sj) for Isajet Yukawa-unified 
solutions with R < 1.05 against uiiq. The experimental central value along with ±2(7 
error bands are also shown. The Wilson co-efficients grow with tan /?, and since we are at 
tan/3 ~ 50 for Yukawa-unified SUSY, we expect large MSSM contributions to BF{b 57). 
In fact, we see from Fig. H that the BF{h — > 57) is typically below the measured error 
bands for mig ~ 10 TeV, and as mig gets very large, the branching fraction approaches 
the SM value. This at least seems to favor mig values ~ 8 — 10 TeV in Yukawa-unified 
models. This observation is also recorded by Altmannshofer et al.[^]. 

In Fig. ^, we show the branching fraction BF{Bs — > evaluated [^] using 



Isatools. The experimental 95% CL upper limit from CDF collaboration is also shown|51|. 
We see that for lower values of mig, in many cases the branching fraction can be near 
or even above the experimental limit. As mig increases, so does m^, which mediates the 
decay. The large value of mie, and hence tua, acts to suppress this branching fraction to 
values about an order of magnitude below present experimental limits. 



3. Axion and axino dark matter 



3.1 Axion cold dark matter 



The axion arises as a by-product of the Peccei-Quinn solution to the strong CP problem|18 



-10- 



zi. 
t 



' '.i •'A.r* it -- 



Si' 





4 6 8 10 12 14 16 

'"is (TeV) 

Figure 6: Predictions for BF{Bs ji^ vs. miQ from Yukawa-unified models generated by 
Isajet7.79 with R < 1.05 and to* = 172.6 GeV. The green dots show results from MCMC scans 
with fixed mig. We also show the measured 95% CL upper limit on this branching fraction from 
CDF measurements. 



2C]. The strong CP problem has its origin in an allowed QCD Lagrangian term 

0^2 



327r2 '^'^ 



(3.1) 



(G^j^is the gluon field strength tensor and G""^^ its dual) which is P and T- violating, but 
G conserving, and hence CP violating. When QCD is coupled to the electroweak theory, 
d is replaced by ^ = + arg{det m^), where niq is the quark mass matrix. The measured 
value of the neutron electric dipole moment (EDM) requires 9 ^ 10^^''. Explaining the 
tininess of this Lagrangian term is the strong CP problem. 

The Peccei-Quinn solution to the strong CP problem invokes a theory with a global 
U{1) (Peccei-Quinn or PQ) symmetry, which is classically valid, but broken spontaneously, 
and by quantum anomalies. A consequence of the broken PQ symmetry is the existence 
of a pseudo-Goldstone boson field: the axion a(x)|jl^. In this case, the Lagrangian also 
contains the terms 

(3.2) 



where we have introduced the PQ breaking scale fa and N is the model-dependent color 
anomaly of order 1. The effective potential for the axion field V{a{x)) has its minimum 
at {a{x)) = —9fa/N, and so the offending GG term essentially vanishes, which solves the 
strong CP problem. A consequence of this very elegant mechanism is that a physical axion 
field should exist, with concommitant particle excitations. 

The axion mass can be computed using current algebra techniques, and is given by 

'""^"^^ (3.3) 



6 eV- 



fa/N 



The axion field couples to gluon-gluon (obvious from Eq. ( |3.2D ) and also to photon-photon 
and fermion-fermion. All the couplings are suppressed by the PQ scale fa- 
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Astrophysical limits from cooling of red giant stars and supernova 1987a require 
fa/N ~ 10^ GeV, or ma ~ 3 x 10~^ eV. In addition, axions can be produced via vari- 
ous mechanisms in the early universe. Since their lifetime (they decay via a — > 77) turns 
out to be longer than the age of the universe, they can be a good candidate for dark 
matter in the universe. Since we will be concerned here with re-heat temperatures of the 
universe Tr ~ 10^ GeV < fa (to avoid overproducing gravitinos in the early universe), 
the axion production mechanism relevant for us here is just one: production via vacuum 
mis-alignment fS^ . In this mechanism, the axion field a{x) can have any value ~ fa at 
temperatures T ^ Aqcd- As the temperature of the universe drops, the potential turns 
on, and the axion field oscillates and settles to its minimum at —0fa/N. The difference in 
axion field before and after potential turn-on corresponds to the vacuum mis-alignment: it 
produces an axion number density 

riait) ~ ^ma{t){aHt)), (3.4) 

where t is the time near the QCD phase transition. Relating the number density to the 
entropy density allows one to determine the axion relic density today: 

An error estimate of the axion relic density from vacuum mis-alignment is plus-or-minus a 
factor of three. Axions produced via vacuum mis-alignment would constititute cold dark 
matter. 

The axion relic density from vacuum mis- alignment, along with error bands, is shown in 
Fig. ^. However, in the event that (a^(i)) is inadvertently small, then much lower values of 
relic density could be allowed (or much higher if (a^(t)) is inadvertently large). Additional 
entropy production at t > tqcD can also lower the axion relic abundance. Taking the value 



of Eq. (3^) literally, and comparing to the WMAP5 measured abundance of CDM in the 



universe, one gets an upper bound fa/N -^^^ 5 x 10^^ GeV, or a lower bound ~ 10 ^ 
eV. If we take the axion relic density a factor of three lower, then the bounds change to 
fa/N ~ 1.2 X 10^2 Qgv, and ~ 4 x 10"^ eV. 

3.2 Warm and cold axino dark matter 

3.2.1 Non-thermally produced axino dark matter 

Since we are working in a supersymmetric model, the axion field will be only one element 
of an axion left chiral scalar superfield 

k = + iV2HaA^) + iOOMx), (3.6) 

where 6 here are the anti-commuting Grassman superspace dimensions arranged in a Ma- 
jorana spinor, and = + ^^757u^[^]. The superfield 4>a contains the i?-even spin-0 
saxion field s, which gets a mass of order the SUSY breaking scale, and the i?-odd spin-^ 
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Figure 7: Axion relic density due to vacuum mis-alignment versus rUa (lower scale) and fa/N 
(upper scale). The plot includes a factor three error estimate along with the WMAP5 CDM 
measured abundance. 



axino field ipa = o,, whose mass is model-dependent, and can range over the keV-GeV 
scale[||, m. 

Here, we assume that the d is the LSP, so that the neutralino is in fact unstable, and 
decays dominantly into x? ~^ ^7- The width T{xi ~^ ha-s been calculated in Ref. | ]23| , 
and is given by 



(1)2 



xi 



2 



,(1) 



128^3 C0S2^T^ Ua/Nf 



where ^4'^'' denotes the bino fraction of neutralino x?) is the model-dependent anomaly 




(3.7) 



factor {e.g. = 1 (6) for KSVZ||5J] (DFSZ||5g]) axions), and CoYY is a model-dependent 
coupling factor [e.g. Cavy = 8/3 in the DFSZ model). 

In Fig. ^, we plot the x? lifetime in seconds versus m^a for four choices of fa/N, and 
taking CaYY = 8/3. The lifetime ranges from lO^^s for fa/N = 10^ GeV, up to ~ 40 
s for fa/N = 10^2 GeV. In the latter case, the will decay while BBN is ongoing. The 
dominant decay into a high energy photon should thermalize with the electron-nucleon 
plasma. It is of note that x? — *■ IQO. three body hadronic decays via intermediate 7 and Z 
can also occur at a small branching fraction. These hadronic decays would be more likely 
to be a threat to disrupt Big Bang Nucleosynthesis. Also, Xi ~^ ^^l can occur, but only 
for m^o > Mz, which rarely occurs in our Yukawa- unified scenario. 

The axino dark matter produced from neutralino decay would compose non-thermally 



produced (NTP) dark matter. Jedamzik et aZ.|26] have calculated the rms velocity profile of 
axino dark matter coming from neutralino decay. A comparison against data from Lyman 
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Figure 8: Lifetime in seconds of the lightest neutralino versus its mass for fa/N — 10^, 10^*^, 10^^ 
and 10^^ GeV, respectively. 

alpha forest leads them to conclude that non-thermally produced axinos will consitute 
warm dark matter for rria ~ 1 GeV. 

The relic abundance of non-thermally produced axinos can be simply obtained from 
the neutralino abundance. The neutralino thermal abundance calculation proceeds by 
solving the Boltzmann equation for neutralinos from freeze-out to the present day, after 
inputting the usual neutralino annihilation and co- annihilation cross sections. Since each 
neutralino decays to one axino, the axinos inherit the neutralino number density, and the 
non-thermally produced axino abundance is simply 

^^TPf^2 ^ I!}^n,h\ (3.8) 

In this regard, if the ratio rria/m^o is small, then large factors of neutralino dark matter 
density can be shed by undergoing Xi ~^ 0,1 decay. It is this mechanism that allows one to 
reconcile the huge neutralino relic abundance from Yukawa-unified models with the WMAP 
measured abundance. 

3.2.2 Thermally produced axino dark matter 

In our scenario, where we only consider Tr ^ /„, the axinos in the early universe are too 
weakly interacting to be in thermal equilibrium. Nevertheless, they can be produced by 
radiation off other particles which are in the thermal bath, much the same as gravitinos 
can be produced in the early universe. Initial calculations of the thermally produced (TP) 
axino abundance were performed in Ref. ||2^ , wherein a variety of QCD axino production 
processes (such as gg — > ag, gq — > dq, • • •) were considered. Divergent diagrams involving 
t-channel exchange of massless gluons were regulated by introducing a "plasmon" mass, 
representing the effective gluon mass in the plasma of the early universe. A later evaluation 
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of thermally produced axino matter in Ref. used the hard thermal loop resummation 
technique of Braaten-Pisarski|5(;] and obtained a reduced axino dark matter yield by a 
factor of ~ 3. The results of Ref. p^] are summarized in the expression 



..TP, 2 rr6, /1.108\ /lO" GeV\ V ma \ f Tr \ , , 

where Qs is the strong coupling evaluated at Q = Tr {e.g. gs = -915 at Q = 10^ GeV from 
our ISAJET RGE calculations). The thermally produced axinos qualify as cold dark matter 
as long as ~ 100 keV|2|, 

In Fig. y, we show bands of the rria vs. Tr plane which give QTP/i2 within the WMAP- 
measured dark matter abundance, for fa/N = 10^^, 10^^ and 10^^ GeV. We see that for 
the lower range of fa/N ~ lO^'' GeV, very low values of ma and Tr are required. In this 
case, with rria ~ 100 keV, the thermally produced axino DM would likely constitute warm 
DM, and furthermore, the low value of Tr excludes some of the possible mechanisms for 
baryogenesis. In this case, if we want Tr 10 GeV and rria ^ 100 keV with dominant 
thermal production of cold axino DM, then we will need higher fa/N ^ 10^^ GeV. 

In Fig. |10|, we show bands of ^J'^h'^ = 0.11, 0.03, 0.01 and 0.001 in the rria vs. Tr 
plane for fa/N = 10^^ GeV. In this case, if the thermally produced axino dark matter 
only constitutes a small fraction of the total dark matter, and most of the remainder is 
composed of cold axions, then much smaller values of rria are allowed, and the thermally 
produced axinos can be either warm or even hot dark matter. 
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Figure 10: Thermally produced axino relic density in the ma vs. Tn plane for fa/N — 10^^ GeV. 



4. The gravitino problem, non-thermal leptogenesis and mixed axion/axino 
dark matter 

4.1 The gravitino problem 

A problem common to all SUSY models including supergravity (SUGRA) is known as 
the gravitino problem. In realistic SUGRA models (those that include the SM as their 
sub- weak-scale effective theory), SUGRA is broken in a hidden sector by the super Higgs 
mechanism. A mass for the gravitino G is induced by SUGRA breaking, which is commonly 
taken to be of order the weak scale. The gravitino mass m^i2 sets the mass scale for all the 



soft breaking terms, so that all SSB terms end up also being of order the weak scale [57, 58 1. 

The coupling of the gravitino to matter is strongly suppressed by the Planck mass, so 
the G in the mass range considered here ~ '^le ~ 5 — 20 TeV) is never in thermal 

equilibrium with the thermal bath in the early universe. Nonetheless, it does get produced 
by scatterings of particles that do partake of thermal equilibrium. Thermal production of 
gravitinos in the early universe has been calculated in Refs. |^^, where the abundance 
is found to depend naturally on and on the re- heat temperature Tr at the end of 

inflation. Once produced, the Gs decay into all varieties of particle-sparticle pairs, but 
with a lifetime that can exceed ~ 1 sec, the time scale where Big Bang nucleosynthesis 
(BBN) begins. The energy injection from G decays is a threat to dis-associate the light 
element nuclei which are created in BBN. Thus, the long-lived Gs can destroy the successful 
predictions of the light element abundances as calculated by nuclear thermodynamics. 

The BBN constraints on gravitino production in the early universe have been calculated 
by several groups |^^. The recent results from Ref. [^] give an upper limit on the re- 
heat temperature as a function of m^^2- The results depend on how long-lived the G is 
(at what stage of BBN the energy is injected), and what its dominant decay modes are. 
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Qualitatively, for 711-3/2 ~ 5 TeV, Tr ~ 10^ GeV is required; if this is violated, then too 
many G are produced in the early universe, which detroy the ^He, ^Li and D abundance 
calculations. For 1713/2 ~ 5 — 50 TeV, the re-heat upper bound is much less: Tr ~ 10^ 
GeV (depending on the ^He abundance) due to overproduction of '^He arising from n ^ p 
conversions. For 7713/2 ^ 50 TeV, there is an upper bound of Tr ~ 10^*^ GeV due to 
overproduction of x? LSPs due to G decays. 

Solutions to the gravitino BBN problem then include: 1. having 7713/2 ~ 50 TeV but 
with an unstable Xi (-R-parity violation and no Tr bound), 2. having a gravitino LSP so 
that G is stable or 3. keep the re-heat temperature below the BBN bounds. We will here 
adopt solution number 3. In the case of 50(10) SUSY GUT models, we expect 7713/2 ~ 77715 
and since 77716 ~ 5 — 20 TeV, this means we need a re-heat temperature Tr ~ 10^ GeV. 



4.2 Non-thermal leptogenesis 

The data gleaned on neutrino masses during the past decade has led credence to a particular 
mechanism of generating the baryon asymmetry of the universe known as leptogenesis [l6l| . 
Leptogenesis requires the presence of heavy right-handed gauge-singlet Major ana neutrino 
states ipNc{^= Ni) with mass Mj^. (where 7 = 1 — 3 is a generation index). The Ni states 
may be produced thermally in the early universe, or perhaps non-thermally, as suggested in 
Ref. via inflaton cj) — > NiNi decay. The Ni may then decay asymmetrically to elements 
of the doublets- for instance r(A'^i — > h^e~) 7^ r(A^i — > h~e~^)- owing to the contribution 
oi CP violating phases in the tree/loop decay interference terms. Focusing on just one 
species of heavy neutrino A^i, the asymmetry is calculated to be|63] 



^^ r(iv.^,.) rw^r) ^ 3M^ ^^^^ 

Fatj Stt 

where 7771^3 heaviest active neutrino, Vu is the up-Higgs vev and dgfr is an effective 

CP-violating phase factor which may be of order 1. The ultimate baryon asymmetry of the 
universe is proportional to e, so larger values of MjVj lead to a higher baryon asymmetry. 

To find the baryon asymmetry, one may first assume that the A^i is thermally produced 
in the early universe, and then solve the Boltzmann equations for the B — L asymmetry. 
The ultimate baryon asymmetry of the universe arises from the lepton asymmetry via 



sphaleron effects. The final answer] 64], compared against the WMAP-measured result 



^ ~ 0.9 X 10 for the baryon-to-entropy ratio, requires M^^ ~ 10^'^ GeV, and thus 
a re-heat temperature Tr ~ lO^'^ GeV. This high a value of reheat temperature is in 



conflict with the upper bound on Tr discussed in Sec. iA_. In this way, it is found that 
generic SUGRA models with i?-parity conservation are apparently in conflict with thermal 
leptogenesis as a means to generate the baryon asymmetry of the universe. 

If one instead looks to non-thermal (NT) leptogenesis, then it is possible to have lower 
reheat temperatures, since the Ni may be generated via inflaton decay. The Boltzmann 



equations for the B — L asymmetry have been solved numerically in Ref. |65]. The B — L 



asymmetry is then converted to a baryon asymmetry via sphaleron effects as usual. The 
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baryon-to-entropy ratio is found to be|65| 



where niif, is the inflaton mass. Comparing calculation with data, a lower bound Tr ~ 10 
GeV may be inferred for viable non-thermal leptogenesis via inflaton decay. 

4.3 A consistent cosmology for Yukawa-unified models with mixed axion/axino 
dark matter 

Next, we investigate various scenarios with mixed axion/axino cold and warm dark matter, 
first to see if they can yield a consistent cosmology, and second, to see if cosmology provides 
any insight into allowed model parameters. Here, we will investigate four cases. 



1. We take fa/N = 10^^ GeV. Using the central value from Eq. 3^, we obtain a small 
fraction of axion CDM: nah^ ~ 0.017. The bulk of CDM must then be composed 
of something else: in our case, thermally produced axinos (so ~ 100 keV). We 
take fiT^/i-^ = 0.083. Then to obtain the WMAP5 measured DM abundance, we get 

oFP/i2 ^ 0.01. 



a 



2. We take fa/N = 4 x 10^^ GeV so that the central value of axion relic abundance Eq. 
U yields Q^/i^ ~ 0.084. This gives dominant axion CDM, so that thermally and 
non-thermally produced axino abundance should be small. Here we assume an equal 
mix of thermal and non-thermal axinos, so Qj^h"^ = Q^'^^h'^ = 0.013. 

3. We take fa/N = 10^^ q^v, and a factor of 1/3 error on ^lah as in the lower dashed 
curve of Fig. 0. Thus, ^ah"^ = 0.084 so we have dominant axion CDM. As in the 
previous case we take an equal mix of thermally and non-thermally produced axinos: 

nTPh^ = n^TP^2 ^ 0.013. 

4. Here, we again take fa/N = 10^^ GeV, but assume the axion vev is accidentally close 
to zero so that it is nearly aligned with the potential minimum, instead of mis-aligned. 
Even though fa/N is large, the resulting axion abundance is small: Qah'^ ~ 0. In 
this case, thermally produced axinos should make up the dominant CDM component. 
We take Qj^h"^ = 0.1 and O^^/i^ ~ 0.01. This case was shown previously as one 
adopted in Ref. ||3^ . 

Once the value of fi^^/i^ is known, we may calculate ^■^oh'^ and m^o in any super- 
symmetric model (with a Xi NLSP) and use Eq. ( |3.8[ ) to calculate the value of ma that 
is needed. Then, if a value of fa/N has been selected, and we know Qj^h"^, we can use 
Eq. ( p. 91 ) to determine the required re-heat temperature Tr. We plot our final results in 
the rria vs. Tr plane, so that we may see whether a consistent cosmological scenario may 
be found for any of our Yukawa-unified solutions, and also whether a consistent cosmology 
helps to select out preferred values of the soft SUSY breaking parameters. 

Our main results are shown in Fig. Here, we generate Yukawa-unified solutions 
with R < 1.05 first with Isajet7.79 through MCMC scans for mig = 5 TeV (small-red 
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Figure 11: Plot of locus in the ma vs. Tr plane of four Yukawa-unified cases of mixed axion/axino 
dark matter, along with four different mie values. The thermally produced axino relic density 



calculation is only valid for Tr 
10^ GeV. 



10 GeV. Non-thermal leptogenesis prefers the region with Tr 



points), 8 TeV (dark-blue larger points), 10 TeV (medium-blue larger points) and 15 TeV 
(very large light-blue points). The cases 1-4 are labelled as C1-C4. For the region below 
Tr ~ 10^ GeV, the calculation of thermally produced axinos breaks down. Moreover, the 
value of Tr is becoming comparable to the sparticle mass scale {i.e. the heavier sparticles 
will not be produced in thermal equilibrium in the early universe) so the calculation of 
0^0 /i'^ would also break down. Values of Tr ~ 10^ are compatible with non-thermal 
leptogenesis via inflaton decay. 

• We see from Fig. ^ that case CI with fa/N = 10^^ GeV and dominantly TP axino 
dark matter gives almost all solutions in the region with Tr < 10^ GeV. The solutions 
with rriiQ ~ 5 — 8 TeV especially have low values of Tr. The solutions with mie ~ 15 
TeV do have Tr ~ 10^ GeV, but these solutions also have < 10~^ GeV, and so 
the dominant DM component from thermally produced axinos is likely warm DM. 
This scenario would thus be difficult to accept cosmologically, for any value of mig. 

• For case C2, we have fa/N = 4 x 10^^ GeV and dominant axion CDM. Here, the 
larger value of fa/N allows solutions with a similar value of rria as case CI, but with 
a higher value of Tr. The solutions with mig = 10 — 15 TeV do emerge with Tr ~ 10^ 
GeV- in the range for NT leptogenesis. Many of these solutions still have nia < 10~^ 
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Figure 12: Plot of locus in the ma vs. Tr plane of four Yukawa-unified cases of mixed axion/axino 
dark matter for all mig values. We show results from ISA jet (blue) and Softsusy (green). 

GeV, so that the thermally produced axinos are warm. In this case, < 10~^ GeV 
is allowed, since instead the axions actually make up the CDM. 

• In case C3, fa/N = 10^2 Qgv with dominant axion CDM. The larger fa/N gets, the 
larger are the calculated values of Tr. The thermal and non-thermal axinos both have 
small contributions to the relic density, so the entire band of solutions with Tr > 10® 
GeV yields a consistent cosmology. This case requires an axino with nia ~ 10^'^ GeV 
for Tr > 10® GeV. 

• Finally, case C4 is constructed to have a large value of fa/N = 10^^ GeV, but with 
a tiny axion relic abundance due to accidental vacuum alignment. In this case, the 
thermally produced axinos comprise the CDM. Solutions are found with Tr > 10® 
GeV for rria ~ 6 x 10~^ GeV. However, in this case, the solutions with rria ~ 10^^ 
GeV would not be allowed, since they likely yield a dominant warm DM scenario, in 
contrast to requirements from large scale structure formation that the bulk of DM 
be cold. 



In Fig. 12, we show Yukawa-unified solutions with R < 1.05 for cases C1-C4 in the 
nia vs. Tr plane from both Isajet7.79 and Softsusy2.0.18 MCMC scans. ^ Isajet results 
are in blue, while SOFTSUSYresults are in green. Both sets of results line on the same line 



'To be precise, the results for mie = 5, 8, and 10 TeV come from both Isajet7.79 and Softsusy 2.0.18. 
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for a given case. The Isajet results reach to somewhat higher Tr values than Softsusy. 
This is due in part because Isajet can more easily generate Yukawa-unified models for 
very high mig ~ 10 TeV, and these models typically have larger ^y^oh"^ values, and hence 
smaller values of rria- Also, the Softsusy results tend to have lower /i values than Isajet 
results, which also tends to lower the value of the neutralino relic density. Thus, case CI 
does not lead to a consistent cosmology for either Isajet or Softsusy. Isajet can obtain 
cosmologies with Tr > 10^ GeV for cases C2-C4, but Softsusy can generate cosmologies 
with Tr > 10^ GeV only for case C3 (just barely) and case C4. 



5. Summary and conclusions 

One vestige of supersymmetric SO (10) grand unified theories may be that the third gen- 
eration t — b — T Yukawa couplings unify, in addition to gauge couplings and matter multi- 
plets. Assuming the MSSM is the low energy effective theory at energy scales Q < Mqut, 
we are able to find parameter space solutions that yield a superparticle mass spectrum 
with Yukawa coupling unification good to 5% or better, using the Isajet 7.79 and Soft- 
susy 2.0.18 programs. The sparticle mass spectrum that qualitatively emerges is that first 
and second generation scalars lie in the multi-TeV regime, third generation scalars, /i and 
VIA lie in the few TeV range, and gauginos lie in the sub- TeV range. The neutralino relic 
density turns out to be 10 — 10^ times the measured dark matter density, prompting the 
suggestion that in this case, the axino is a better LSP candidate, so that the dark matter 
of the universe would be composed of an axion/axino mix. 

Yukawa-unified SUSY models should thus give rise to three components to the dark 
matter density: axion dark matter produced via vacuum mis-alignment at the QCD phase 
transition, non-thermally produced axinos from Xi ~^ ^7 decay (likely warm dark matter 
if rria ~ 1 GeV), and thermally produced axinos which are likely cold dark matter unless 
nia ~ 100 keV. We compute the abundance of all three components of dark matter in 
Yukawa-unified models for four different scenarios containing either dominant axion or 
dominant axino cold dark matter. The relative abundances depend on J7^o/i^, rn^o, rria, 

fa/N and Tr. It is found that for all solutions with m-ig ~ "1-3/2 ~ 5 TeV, the value of Tr 
needed is below limits calculated from BBN constraints, thus solving the gravitino BBN 
problem. 

We also find that it is very difficult for models with PQ symmetry breaking scale 
fa lower than ~ 2 x 10^^ GeV (case CI) to generate dominantly cold dark matter and a 
sufficiently large Tr to be consistent with at least non-thermal leptogenesis, for any allowed 
value of > 5 TeV. However, if the PQ scale is large enough (~ 4 X 10" GeV), then 
we can generate a universe with axions as the dominant component of CDM (cases C2 and 
C3), and only a smaller component is composed of possibly warm axinos with rria ~ 10^^ 
GeV. These solutions work out if mig ~ 10 TeV, so that the higher side of the mie ~ 3 — 15 



TeV range is preferred. This is also in accord with results presented in Sec. 2.3, where it 



Points with mie = 15 TeV come only from Isajet 7.79, because Softsusy 2.0.18 does not give consistent 
EWSB for such high mie- 
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is found that BF{b —>■ 57) and BF{Bs — > ^"'"^~) decay also prefer mig ^ 10 TeV. Finally, 
we show a case C4 with accidentally small axion CDM component, but which gives cold 
thermally produced axino dark matter, as long as rria ~ 100 keV. This scenario allows for 
solutions with mig as low as 5 TeV. 

As far as tests go of the Yukawa-unified SUSY models with mixed axion/axino dark 
matter, it is clear from previous work that very large signal rates from gluino pair pro- 
duction and subsequent cascade decays should be visible at LHC relatively soon after 
start-up [|^]. The x? produced in cascade decay events will still yield events with i?™*^^, 
since the — ^ ^7 decay occurs far outside the detector. In cases C2 and C3, there is a good 
chance for direct detection of relic axion dark matter at experiments such as ADMX [l67| . 
However, direct and indirect searches for WIMP dark matter would likely turn up null 
results. 



A. Appendix: Yukawa-unified benchmark points from Isajet7.79 

Here, we present several updated benchmark points from Isajet7.79 with Yukawa-unified 
solutions. Points A, B and C with mig = 5, 10 and 15 TeV respectively all require mixed 
axion/axino dark matter. Point H has mie ~ 3 TeV, and can accomodate neutralino dark 
matter since neutralino annihilation through the h resonance leads to a neutralino dark 
matter abundance in near accord with WMAP5 measurements.^ 
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